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Abstract 

The one-particle inclusive cross section in deeply inelastic lepton-nucleon scattering, expressed in 
terms of parton densities and fragmentation functions being differential in the invariant mass of 
the observed hadron and of the incoming nucleon, diverges if this invariant mass vanishes. This 
divergence can be traced back to the kinematical configuration where the parent parton of the 
observed hadron is emitted collinearly from the incoming parton of the QCD subprocess. By 
using the concept of "fracture functions" , which has recently been introduced by Trentadue and 
Veneziano, it is possible to absorb this divergence in these new distribution functions as long as 
the observed hadron is not soft. This procedure allows the determination of a finite one-particle 
inclusive cross section in ncxt-to-leading order QCD perturbation theory. We give details of the 
calculation and the explicit form of the bare fracture functions in terms of the renormalized ones. 
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1 Introduction 



QCD as a theory of the strong interactions should be able to give a description of one-particle inclusive 
processes. Since non-perturbative effects such as the nucleon wave functions and the fragmentation of 
quarks and gluons are not manageable in explicit calculations, it is necessary to parametrize them by 
means of phenomenological distribution functions in order to make contact with perturbative QCD 
[p], H|, ||, |4], ||, ||, The case of one-particle inclusive lepton-nucleon scattering has been considered 
by several authors ||, |9], S, ||, ITO, 11 , Jl2|, |DJ . The canonical procedure to do the calculation is to 



define a variable z^ = P Q h/P q as in H, where P is the momentum of the incoming nucleon, q is 
the momentum of the exchanged virtual photon, and h is the momentum of the observed hadron. 
It turns out that da/dz^, the differential cross section for the production of h plus anything else, 
is singular for — > 0. A similar effect (da/dz^ — » oo for z^ — > 0) can be observed in the case of 
e + e~ annihilation, where z^ —> means that the observed hadron h is soft (E^ — > 0). This singularity 
is physical, being related to the fact that the multiplicity for the emission of soft particles is infinite. 
However, in lepton-nucleon scattering the situation is different: z^ — > means (in the case of massless 
particles) that either E^ — > or that h is collinear with the target remnant. In the case E^ — > the 
same argument as in the e + e~case applies. If, however, z^ — » 0, but not Eh — > 0, the process under 
consideration is that of a hadron h with non-vanishing energy emitted in the remnant direction, which 
is physically perfectly well defined and has to be distinguished from the case of a soft process. It is 
possible to separate the two cases by introducing energy and angle variables z and v {z = means 
£"h = 0, v = 1 means that h is collinear with the remnant) for the observed hadron. This allows to 
consider the case of z^ — > without having E^ — > 0. da/dzdv can be calculated in perturbative QCD 
by means of parton densities and fragmentation functions. Again, it turns out that da/dzdv — » oo for 
z 7^ 0, v — > 1. A closer analysis of the problem reveals that it stems from the contributions in which 
the fragmentation function is attached to an outgoing parton which is emitted from the incoming 
parton. It is clear that a collinear divergence appears if the mentioned outgoing parton is collinear to 
the incoming one. In this phase-space region, the observed hadron h comes from the target remnant 
region. This singularity cannot be absorbed into the parton distribution functions or the fragmentation 
functions since in lowest order (in the formalism based solely on parton densities and fragmentation 
functions) the observed hadron always originates in the current region. To absorb the divergence, it 
is necessary to introduce a distribution function that describes hadrons in the target remnant region. 
This has been done in [O]; the new distribution functions have been called "fracture functions". In 
this paper it is shown that, using these concepts, a finite cross section da/dzdv for deeply inelastic 
lepton-nucleon scattering can be calculated for E^ ^ in next-to-leading order QCD perturbation 
theory. 

The outline of this paper is as follows. In Section [2| the calculation of one-particle inclusive cross 
sections is briefly reviewed. The 0(a s ) QCD corrections to the parton-model process are determined 
in Section [3|. In Section |I] the explicit form of the renormalized fracture function in terms of bare 
distribution functions is stated. The finite cross section is discussed in Section ||. The paper closes with 
a summary and some conclusions. Technical details and explicit formulas are relegated to appendices. 
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2 One-Particle Inclusive Cross Sections 



In the following the scattering process 

l + A-fl' + h + X (1) 

is considered, where 1 and A are the incoming charged lepton and nucleon, respectively, 1' is the 
scattered charged lepton, h is the identified hadron, and X denotes anything else in the hadronic 
final state. Let I and I' be the momenta of the incoming and outgoing lepton, respectively, P the 
momentum of the incoming nucleon and h the momentum of the outgoing identified hadron. The 
integration over the angles that determine the relative orientation of the leptonic and hadronic final 
state can be performed. The remaining lepton phase-space variables are the Bjorken variables 

Q 2 PoQ 

Here q = I — I' is the momentum transfer and Q 2 = —q 2 > is the virtuality of the exchanged virtual 
photon. 

The cross section for the production of an n-parton final state differential in x B , y and the phase-space 
of the final-state partons is 

- E/f W (y"W + y'% (3) 

Here the case of one-photon exchange has been considered; the exchange of weak vector bosons is 
neglected for simplicity. S H = (P + I) 2 is the square of the total CM energy of the lepton-nucleon 
scattering process, £ is the momentum fraction of the initial parton of the QCD subprocess, and 
Pi/a(0 is the probability distribution function for parton i in the nucleon A; Pj/a(£) can be either a 
parton density /j/a(0 or a fracture function Mj/j/a(£, C)j see below, a = e 2 /Air is the fine structure 
constant and d = 4 — 2e is the space-time dimension used in dimensional regularization. A factor of 
1/4 for the average over the spin degrees of freedom of the incoming particles is already included. The 
last factor, H^, is the hadron tensor defined by 

= M^M„, (4) 

spins 

where e M (A)A4 M is the matrix element for the process 

7* + parton — > n partons, (5) 
with e M (A) the polarization vector of a virtual photon with polarization A. The ratios 

r"= 1 + V> 2 and yL = 4(1 -») + ! + (!-»)' (6) 

specify the y-dependence of the contributions from the two photon polarizations. The projections 
operating on are the result of the integration over the angles that describe the lepton-hadron 
orientation. The cross section consists of two parts, proportional to Y M , the "metric" contribution 
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(it is obtained by a contraction of the hadron tensor with (—g^ u )), and proportional to Y L , the 
longitudinal contribution. 

For the production of an identified hadron h in the final state, there are in principle two mechanisms 
[ffq]: h may originate from the current (in the parton-model, from the produced quark, or in the 
QCD-improved parton-model, from one of the produced partons of the QCD subprocess) or from the 
target remnant 0], so that 

CT — "current G target- (7) 

These two contributions cannot be disentangled in the QCD-improved parton-model, because c curren t 
gives a contribution to cr targe t in the case of collinear singularities in the initial state. This fact is 
reflected in the definition of the renormalized fracture functions, which are, contrary to the case of 
parton densities, inhomogeneous as a function of the bare fracture functions, as can be seen in Section ||. 
In the QCD-improved parton-model, cr curren t and crtarget are given by the graphs in Figs. 1 (a) and 
1 (b), respectively; cr C urrent can be described by means of parton densities /j/a(0 and fragmentation 
functions D h j i {z). The lowest-order contribution is given by the graph in Fig. 2 (a), crtarget) however, 
needs a new phenomenological input, the "fracture function" M^m^, C)> giving the probability to 
find a parton i with momentum fraction £ and a hadron h with momentum fraction £ in the nucleon 
A, for the parton model process see Fig. 2 (b). The kinematical restriction on £ and £ is £ + £ < 1. 

It is reasonable to consider the production process in the CM frame of the incoming nucleon and 
the incoming virtual photon, so P + q = 0. The positive z-axis is defined by the cy-direction. The 
hadron h has polar angle i? relative to the virtual photon and energy in this frame. The energy of 
the incoming nucleon is 

= % , ,! r (8) 

2 y/X B {l ~ Xg) 

Two new variables v, z can be defined byEI 

The range of these variables is v, z € [0, 1] if the masses of all particles are neglected. The differential 
cross section calculated in the following is 

d<r(l + A — > 1' + h + X) 

dx B dy dzdv 

It turns out that QCD-corrections to the lowest-order process require subtractions in the collinear 
region that make this differential cross section a distribution instead of a function of the variable v. 
It is therefore more reasonable to consider an observable A(v) and to integrate over v 

(A) = £dv^A{v). (11) 



1 The variable z has no relation to the variable % mentioned in the introduction. 
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This way one defines the expectation value 



A 



dx B dy dz 



(12) 



Explicitly, it is given by 

j ■ 



2S H x B e 2 (2ir) d 



Y M {-gn+Y L ^r 



H 



' N 
%A(f|)(l-*W) + / f /A(i)E^ (| 



U ^ ^(l- x B )A(vA. (13) 



The notation in this formula is the following: N is the number of final-state partons in the QCD 
subprocess, u = x B /£, i is a multi-index specifying the parton configuration, so i = (h,i2, ■ ■ ■ , fw)> 
where i a € {u,u, d,d, . . . ,g}, p = (pi,P2, ■ ■ ■ ,Pn) is the set of momenta of the outgoing partons, E a 
is the energy of the a th parton in the (P + q = 0) frame, and v a = (1 — cos$ a )/2, where d a is the 
polar angle of the a th parton in the same frame; M, f and D are the fracture functions, the parton 
densities and the fragmentation functions, respectively. 



It is convenient to define the following constants and projection operators: 



or 



2S B Xg 



2tt-4(1-€) Ql 



L — 



4x B 2 

Q 2 



pp.pv 



(14) 
(15) 



Here, eQi is the electric charge of the quark of flavour i. 

The leading order given by the processes of the parton- model depicted in Figs. 2 (a) and 2 (b) is 



-4lo 



yM ^ 



=q,q 



Wi,h/A ( V' (1 ~ Xb)z ) 5(1 " u) (1 " Xb) A( - l) 

+ £t/7^(?)^(^)^-^(i-,)W 



i-(i-«b)* 



(16) 



To this order, there are no longitudinal contributions. They arise in the QCD-improved parton-model 
in 0(a s ), which is considered in the next section. 



3 The Cross Section to 0{a s ) 

This section gives the details of the calculation of the 0(a s ) corrections to the parton-model processes. 
The virtual 1-loop corrections to the leading-order QCD subprocess are shown in Fig. 3, and the real 
corrections in Figs. 4 (a),(b). 
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The overall effect of the virtual corrections is to multiply the leading-order cross section by [16] 



1 + 



a s /47r/x 2 V T(l - e) 



so 



2tt \ Q 2 J T(l - 2e) 

-4.LO+virt. = y M ^2 Ci 



c F _ 2 _L_ 3 i-8-^. 



(17) 



i=q,q 
1 



— M iih/A (^,(l , r ,):) ,S(l-„)(l -, r; ).l(l! 



+ 



i-(i-^s) z 
1 du /"dp 
u J p 



■I'D 



fi/A 



D 



1 + 



u 

rri r 



2tt I Q 2 / T(l - 2e) 



h/i 



<5(1 - «) 8(1 - p) A(0) 



1 1 

-2-^-3- 



IT 

y 



(18) 



where p is the renormalization scale, a s = a s (p 2 ), and C F = 4/3 is one of the Casimir invariants of 
the colour gauge group SU(3). 

The double and single poles in e represent an infrared divergence, which is cancelled by a contribution 
similar to the real corrections but of opposite sign. 

The projections of the hadron tensor for the real corrections are (see, for example, [jl 
1 



PtfH^tfq - qg) 



e 2(2TT) 2d M " r 
.2e 



8TT-^p Ze 2irC F Q z q 4(l-e) 
Ztt * 



(l-e)P + ^ + 



2Q 2 s % , 



+ 2e 



ig°qg 



(19) 



1 



e 2 (2vr) 



2d 



= 8n2Lv*2nC F QU(l-e) 

ZTT 

e 2^2d P M R ^9 qq) 



u 2 1 



4 Q2 2 S «? 



(20) 



a 



87r 2^ M 27r 2 Q ^ 4(1 " e) 



2vr 2^ 9 



1 2Q 2 s 



99 



2q °%q 1 ^ SiqSiq 



1 - e 



(21) 



e 2 (27r) s 

= 8rr JV e 271-^4(1 -e) 



Q 2 



(22) 



The invariants are = 2papb, and p g , p^, p 9 are the momenta of the incoming parton (quark or 
gluon), an outgoing quark, an outgoing antiquark and an outgoing gluon, respectively. The formulas 
already contain the appropriate factors for the average over the colour degrees of freedom of the 
incoming partons. An additional factor of 1/(1 — e) has been provided for the terms with an incoming 
gluon, because gluons have 2(1 — e) helicity states in 4 — 2e space-time dimensions compared to only 
two in the case of quarks. In order to perform the phase-space integrations, suitable parametrizations 



dPS (2) = (2vr) d n 



I) 

d-l 



(23) 



of the two-particle phase-space are needed. The integration variable has been chosen such that it is 
the one that is actually used in the factorization of the collinear singularities. The integration over 
irrelevant angles can be performed. The results for three phase-space parametrizations A, B and C in 
terms of the variables p, w, v! are written out in Appendix [A]. 

By means of these phase-space parametrizations, the integrations can be performed and the infrared 
and collinear singularities can be factorized. The details of this straightforward but cumbersome cal- 
culation are not given here. The results are presented in a form that explicitly shows the cancellation 
of the collinear singularities by means of a renormalization of the phenomenological distribution func- 
tions. The contributions are denoted by Bf , Bf , C M , B{ , B% , C L . Here the indices in B A and C A 
stand for 



a: phase-space region for the integration variable u, 



a = 1: u £ 
a = 2: u £ 



! ' x B + {l-x B )z 



,1 



xb+(1-x b )z 



A: polarization of the virtual photon, 

A = M: metric contribution, 

A = L: longitudinal contribution. 



The observed hadron in the terms B A originates from one of the outgoing partons via a fragmentation 
function, the one in the terms C A comes from the incoming hadron via a fracture function. The B A 
and C A terms receive contributions from the graphs from Figs. 4 (a) and (b), respectively. Throughout 
the calculation the MS factorization scheme both for the parton densities and fracture functions and 
for the fragmentation functions is used. The choice of the factorization scheme defines the finite parts 
unambiguously. One has to choose two factorization scales, one (M"j) for the renormalized distribution 
functions / and M for partons in the incoming nucleon, and one (Mp) for the fragmentation functions 
D. 

The explicit results for the cross sections are collected in Appendix |B|. One sees immediately that the 
IR singularities proportional to 2/e 2 + 3/e cancel among virtual and real corrections. The cancellation 
of the remaining collinear singularities is discussed in the next section, and the finite terms from 
appendix [l| are discussed in Section ||. 



4 Renormalized Distribution Functions 



The renormalization of the bare parton densities /^m and bare fragmentation functions D^/i is done 
in the MS factorization scheme in the standard way ]i~6| , fif, as required by the total cross section 
in lepton-nucleon scattering and the one-particle inclusive cross section in e + e~ scattering. The 
expressions for the bare densities in terms of the renormalized ones for the factorization scales Mj 
and Mp are 



fi/A(0 



Dh/i(t) 



du 



du 



Sij 5(1 - u) + 



la, T(l-e) f Att/i 2 
72^r(l - 2ej 



Sij 8(1 



u) + - 



r(i 



e 2vrr(l - 2e) 



M) 

' ilTfl 2 



Pi 



f 



i/A 



Mi 



D 



~,M 2 D 
u 



(24) 



(25) 



The expression for Mj h/A(£> C) m terms of the renormalized quantity Af^M^, Cj Mj) can be obtained 
by the requirement that all collinear divergences that are not already absorbed into f r and D r are 
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absorbed into M r . By an appropriate definition of the finite parts, the bare fracture function in terms 
of the renormalized fracture function in the MS factorization scheme is 

,2 s 




+ / e+C -, 

- u x. 



(26) 



We have implicitly assumed that repeated indices are summed over. These expressions are valid to 
0{a s ). The formula for M is inhomogeneous. The homogeneous contribution is just the standard 
Altarelli-Parisi evolution caused by the emission of collinear partons from the parton emanating from 
a fracture function. The inhomogeneous term arises because the observed hadron may originate from 
a fragmentation function of an outgoing parton collinear to the incoming one. 

By requiring that the bare fracture functions do not depend on the factorization scale and consequently 
d 



dbiAfJ 



M i)h/A (£,C) = 0, 



(27) 



the evolution equation for -MJ\y A (£, £, M| ) can be derived. It coincides with that in 

Rewriting the bare distribution functions in terms of the renormalized ones in all cross-section formulas, 
and keeping only terms up to 0(a s ), one can see that all collinear singularities cancel. The result is 
a finite cross section to 0(a s ), which is discussed in the next section. 



5 Finite Contributions 

After the cancellation of the infrared and collinear singularities the resulting cross section is finite. 
Because of the subtractions, it is of the form of a convolution of a distribution with phenomenological 
distribution functions. The "+" prescriptions used are defined by ]l8fl: 



b 

dx D +x[ a, b] (x) cp(x) = I dxD(x)((p(x) -ip(a)). (28) 



The variable and the range of the integration are indicated as a subscript; furthermore, the subtraction 
point is underlined. These distributions are used in the expressions of the ^^/ M , which are presented 
in Appendix [(]. The expectation value of the observable A is, differential in x B , y and z: 



dx B dydz ~ "^LO+virt. + 4;' + + ^C M + + + ^C L > ( 29 ) 

with 

A f - Y M V c- 

-^LO+virt. — 1 / t c i 

i=q,q 



{ t m T M WA (v' (1 - Xfl)z ' M /) S(l-n)(l-x B )A(l) 

K \-(l-x B )z 

+ £ t / 7 (t ■ M >) (? M °) k-w-,) m 



a, ( „ 7r 2 



1 + i^-T <30) 
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A f 



V M ° s f xb+(i-*b)z du f 1 dp 



Xg «■ J a(u) P 

f ^(l7' M f) D ^(^ M 'd) \-^P q ^q(u)8(l-p) + C F <S> 



i=q,q 

' Xb -,m 2 

u 



M 
lqq 



M 2 



+/^ A ^,M / 2 ) Dl /g ^Ml) \-\n-^P gq ^ q (u)5(p-a(u)) + C F ^ qg 



Xb -,m 2 

u 

m 



i 



In ^ Pq q ^ g (u) 5(p - a{u)) - In P q ^ g (u) 5(1 - p) + ± $>% q 



(31) 



A 



i=q,q 



f>4^ M i) D 'v(? M °) 



Mf M 2 
In P q ^ q (u) 5(1 - p) - In P q ^ q (p) 5(1 -u) + C F >5> 

Ml 



M 
2qq 



Y M y j c, l ^[ 1 ^A(l) 

«=<?,<? l-(l-a: s ) Z 

+Kh/A ( f , (1 " a*)*, M?) { " ^ ^ W«) (1 - x fl ) + 1 *j 



(32) 



(33) 



A 



Ea s f x B +{i-x B )z du f 1 dp 
°i7r — —A(v(p,u)) 



i=q,q 



L 

lqq 



(34) 



B 2 



= y 



^ 27T 7 -5H U J z p ^ VP ' " 
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L 

2qq 



+fI/A[-,M) 



D 



P 



-,M 2 D ) C F $k 



2<?9 



•I'd 



2 <1> -"" 



(35) 



Al 2vr J m u v ' 



i-(i-^b) 2 



2 9 



(36) 



For some of the contributions the integration region is split into two parts, according to finite terms 
that arise from singular regions. The form given here should be suitable for numerical evaluation. 



6 Summary and Conclusions 

In this paper the one-particle inclusive cross section in deeply inelastic lepton-nucleon scattering has 
been expressed in terms of parton densities, fragmentation functions and fracture functions. A cut in 
the transverse momentum of the observed particle is not required. Compared to the case of a hadron 
that is well separated from the target remnant, additional collinear singularities are encountered 
because of the enlarged phase-space available to the parent parton of the observed particle. They 
stem from the phase-space region where the parent parton of the observed particle is collinear to 
the incoming parton. These singularities can, to 0(a s ), be consistently absorbed into the fracture 
functions. A check of the calculation is provided by the determination of the evolution equation of 



the fracture functions; it coincides with the one given in 14]. 



Fracture functions allow for a coherent and unified description of one-particle inclusive production 
processes in the case of incoming hadrons, without having to assume specific models for the target 
fragmentation fll9f . As in the case of parton densities, the ignorance of the nucleon wave function 
is parametrized in a phenomenological distribution function, which has to measured. For target 
fragmentation, fracture functions have the advantage over phenomenological models that the correct 
scale dependence is built in, because of their evolution equation. It is expected that fracture functions 
are as universal as parton densities, i.e. that they can be determined in a specific process and then 
used for predictions in other processes. 

The results of this paper should be consistent to 0{a s ) if the fracture functions are extracted at a scale 
Mg and then used for predictions at a scale M 2 , with In (M 2 /Mq) not too large. However, in order to 
evolve the fracture functions to scales M 2 very different from Mq , the evolution kernel (the analogon 
of the anomalous dimensions in the case of parton densities and fragmentation functions) should be 
known on the two-loop level. This effort would be necessary to have genuine 0(a s ) predictions in this 
case. 

An interesting question would be to find out whether factorization of QCD subprocesses is always 
possible that collinear singularities can consistently be absorbed order by order into renormalized 
distribution functions, in the case that there is no transverse momentum cut for the observed particle. 
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This would certainly be necessary to make the cross section well-defined. Another question is whether 



there is a connection with the formalism of Mueller [[20| , 21, 22] for one-particle inclusive processes. 
We hope to return to these questions in the future. 



A more practical issue is the application of the concepts of |l3 and the results of this paper to a real 
experimental situation, e.g. to deeply inelastic electron-proton scattering at HERA. An interesting 
application would be the "heavy-quark inclusive" cross section, i.e. the cross section for the production 
of a heavy quark plus anything. There is no obvious reason why the methods used in this paper should 
not be applicable to the production of heavy quarks instead of hadrons. In fact, explicit results for 
the fragmentation function of partons into heavy quarks are available [23], and similar results should 
be possible for fracture functions. Work along these lines is in progress. 
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A Phase-Space Parametrizations 



In this appendix three different parametrizations of the two-particle phase-space are given. 

A) The integration variable p is the energy E\ of one of the partons in the (P + q = 0) frame, scaled 
by the proton momentum, 

L (37) 



Eo(l-x B ) 



This parametrization is used for those contributions in which a collinear singularity has to be absorbed 
into a fragmentation function D; it is 



dps (2) _ 1 W , 2 y* u(l-x B ) 
d b " 8vr IY1 -e)\y) u-x„ 



8vr T(l - e) ) u - x B 
(1 - x B )~ 2e u~ e (1 - u)- e (u - x B ) 2e (p - a{u))~ € (1 - p)~ e dp. (38) 



Here 



a(u) = ^L^. (39) 
1 — x B u 

The energies and invariants are given by 

E 1 = E (l-x B )p (40) 

E 2 = E {l-x B )(l- p + a(u)) (41) 

su = Q 2 — i 12) 

sa = Q 2 ^-^(p-a(u)) (43) 



1 - 


u 


u 




1 - 


X B 


u — 


X B 


1 - 


x B 


u — 


X B 



Si2 = Q 2 ^-^(l-p). (44) 
u — X B 

The angular variable v\ = (1 — cosi?i)/2 is 

x s (l-u)l-p . . 

ui = u(p,«) = . (45) 

u — x B p 

The range of integration is restricted to p G [a(u), 1]; however, in order to ensure that the energy of 
the parent parton of the observed hadron is larger than that of the hadron, the additional condition 
p > z must be satisfied. 

B) The integration variable w is related to an angular variable in the CM system of the virtual photon 
and the incoming parton. Its relation to p is 

w = -}—?- (46) 
1 — a(u) 

This parametrization is used for the contributions that involve a fracture function; it is given by 

dps(2) = hw^T) ( Q2 y e ( i -«)" e « e K i -«')r ed «'- ( 4r ) 

ii 



The energies and invariants are 

E 1 = E (l - x B )(l - (1 - a(u))w) 
E 2 = E (l -x B )(a(u) + (1 -a(u))w) 

S\2 = Q 

U 

sn = Q 2 -(l-w) 

u 

Si2 = Q 2 ~W. 

U 

The angular variable v\ is 

a{u)w 



v\ = v(w, u) = 



1 — (1 — a{u))w 
The range of integration is restricted tow€ [0, 1]. 

C) This parametrization is convenient for the contributions in which the observed hadron orig 
from a parton that is collinear to the incoming parton. A variable v! is introduced by 

/ -, 1 — X B 
u = 1 pu. 

Xg 

The parametrization is given by 

dPS (2) = j_ f Q 2\- e _gg_ 

& 8vr r(l - e) ) u -x B 

■ (u - u')- £ (1 - u)~ £ ( 1 XB 1 ~ U ' ) (u-x B ) 2e x B - e (l-x B )- £ du'. 

\ 1 - x B u J 

The energies and invariants are 
1-u' 



E x = E x B 

u 

E 2 = E {l-x B )(l 
S12 = Q 



x B u — u 
1 — x B u 



u 



o X B u U 

sn = Q 



u — x B u 



2 1 X B I \ X B 1 u ' 
Si2 = Q 1 - ~, 

u — x B \ 1 — x B u 
The angular variable v\ is 



, 1-u (1 - x B )u - x B {l - v!) 

vi = v(u , u) - 



u — x B 1 — u' 

The range of integration is restricted to v! G [1 — u, u]. It has to be further restricted by 
u < 1 z u 

X B 

in order to avoid that the outgoing hadron has a larger energy than its parent parton. 
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B Summary of Cross-Section Formulas 



In this appendix the expressions for the cross section including the divergent parts are written out. 
They are 

A Bf = Y M £ 

i=q,q 



[ f *b+(i-*b)z dn f 1 dp 

\ Jx B U Ja(u) P 

fi/A ) D h/i [~ p 



A(v(p,u)) 



1 T(l -e) ( iirp 2 



P q ^ q (u)5(l-p) 



e T(l - 2e) \ M| 
- In *(l-p) + C F ^j 



In P^,(u) <5(p - a(u)) + C F <&g fl 



| lT(l-e) / 4vr^ 2 

- In - J- jV_ s («) <5(n - «(»)) - In W«) S(l - p) + - *g,J 



r ( i- 2£)V #J *W«W-'> 



+ 



2£j + (l — xg)z ClU 



■I'D 



fi/A 



u 



(l-s B )A(l) 



D 



h/9 



(I-iJzuW 1\ r(l-e) ( 4irn 



1 u 



x B (l-ti)/V T(l - 2e) ^ M| y 1 - u x B 



.-xb\ n ( {i-x B )zu \ r(i-e) /yy 

+/ 9 /A ( u j Mi/i ^ Xs(1 _ n) J ^ J r(l - 2e) V Mj ; 1 - « x s ' ""• " 



+0(e) 



(63) 



A B M 



rM 



i=q,q 



2vr 



B. U Jz P 



x b+( 1 ~ x b) z 
z 



h/i 



i r(i - 6) 

eT(l-2e) 
M 2 



AirpT 



P q ^ q (u)5(l-p) + 



4irp 2 



P q ^ q (p)5(l-u) 



In P q ^ q (u) 5(1 - p) - In ^§> P q ^ q (p) 5(1 -u) + C F J 
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+fi/A [ f ) D h/g (J 



1 T(l -e) 4np 2 



e T(l - 2e) ^ Mj 

Ml. 



P g ^ q (p)S(l-u) 



-ln-^-P g ^ q (p)5(l-u) + C F ^ qg } 



+fg/A [ ) D h/ , 



i „i 

-m-^p q ^ g (u)6(i- P ) + -^ 



+0(e) 



= y m Q — 

A_ 2- .. 

*=<?,<? i-(i-^s)^ 

M hh/A (^,(l-x B )z^ 

- In ^ P^,(u) (1 - x fl ) + C F *Mj 



-I^(^)'w-)d-*) 

1 »,rl 

-ln-^P^ 5 H(l-x fl ) + -$f| 



+0(e) 



ct s f x B +(i-x B )z du f 1 dp 



i=q,q 



U Ja(u) P 



A(v(p,u)) 



L 

19S 



199 



+ 0(e) 
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du f 1 dp 



l=q,q 



+/va ( v 1 /} 



h/i 



p 

z 

p 

z\ 1 



U J z p 



A(v(p,u)) 



C F $ L 



2qq 



+f g/ A -r MV* - 



p7 2 " 2ot 



29.9 



(67) 



Y L J2 Ci— f 1 ^ 

1=9,9 



A(l) 



M, 



i,A/A 



■r-r 



II 



,(l-x B )z^j C F <S> 



+M gA/A [—,{l-x B )z) -$ 



■it 



2 5 



+ 0{e) 



Here the subtracted and unsubtr acted Altarelli-Parisi splitting functions are [16, 24 1: 



P,^,(u) = C F 



Pqq^giu) 



1 

2 

1 



1 - 2u + 2u 2 

2- - 2 + it 
it 



1 



1 - u 
l-2u + 2u 2 



1-u 



(68) 

(69) 

(70) 
(71) 
(72) 
(73) 
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C Explicit Formulas for the Finite Contributions 



In this appendix the explicit expressions for the finite contributions are presented. They are 

.J 



M 
lqq 



7T 



6(1 -u) 5(1- p) — 



+5(1 - p) 

+5(1 - u) 
1 



+t*[0,l] 



+2 



ln(l - u) 
1-u 

ln(l - p) 



1-P /+p[0,l] 
1 



,-, \ i /-, \ 1 + u 2 u — x B 
+ 1 -u - (1 + u)ln(l - u) - In 



1 — u 1 — x B 



1 + P 2 

+ l-p-(l + p) ln(l -p) + - T -!— In p 



1-P/+ P [0,i] VI 



u 



-«[0,1] 



+«[0,1] 



1-P 

(l + ^-f^--) (l + «) 

V 1 - P/+p[o,i] 



+( l- p) ^^fl + ^ 1 -^)V 2 ^^ + 2 



u — X B 



u — X B 



u — X B 



1<29 



= <5(p-a(u)) 



M 



1 + u 2 1-u 
1 — u + m 



1-u 



7/ 



+ 



1 + u 2 (1 - x B ) 2 u 2 x s (l -x B ) l~x B u 2 
+ ~ —- P ; —u - 2- 



p - a(u) J +p [ ( tt ) j i] 1-U (u- X B ) 2 1-U (u - X B ) 2 " II - ;v ! - „ 

1 - ul 



5(p-a(u)) 



2u(l -u) + (l-2u + 2u l ) In ■ 



u 



+5(1 - p) 
1 



2u(l-u) + (l-2u + 2u 2 )ln 



(1 - 2u + 2u 2 ) + 



2m_ (1-u)(1-x b ) 



u — x B 



P - a(u) J + „[ (u),i] 



1 -P/+ P [o,i] 



(1 - 2u + 2u 2 ) - 2 



1 - x B 
u — x B 



-u 



M 
2gg 



M 
1?? 



M 
2gs 



+ 



5(1 -u) 
1 



1-u 



p+ (lnp + ln(l - p)) [p+--2 

(p+--2 

+u[o,i] V P 



+2 - 2^^- + -^p - Xb } 1 * B l U (l -p)+ ' 



1 + u 2 



1i 37r U Xf> 



(u - X B ) 2 



2 1. 1 



p — a(u) 1 — u pi — u 



M 
2<?g 



5(1 -p) 
1 

p — a(u) 



2n (l-u) + (l-2n + 2u 2 )ln (1 " )(J '' d) 



u — x B 



(1 - 2u + 2u 2 ) + 



1 



1 -P/+ P [o,i] 



(1 - 2u + 2n 2 ) - 2 



1 - x B 

u — x B 



-u 



(1 - X B ) 



Z«5(l-u) -ij f _L_ +2 
2 2 U-«/+«[oa] 



ln(l - u) 

1-U / +u [ ,l] 



+3-u-(l + u) ln(l-u) 



1 + u 2 
1-u 
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Inn 



$f = (l-x B ) (hil—H-l) (l-2u + 2u 2 ) (81) 

^ = 4 ^- 2 (!--) (84) 

^ = *f M (85) 

*2 W = (86) 

*2 OT = (87) 

$J = (1 - a«) « (88) 

$J = (1 - x B ) 2«(1 - u) (89) 
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Figure Captions: 

Fig. 1 Contributions to the one-particle inclusive cross section: current (a), target (b). / stands for 
the parton densities, M for fracture functions and P for any hard process. 

Fig. 2 The parton-model process for the current contribution (a) and target contribution (b). D is 
the fragmentation function. 

Fig. 3 The virtual corrections to 0(a s ). 

Fig. 4 The real corrections to 0(a s ): diagrams contributing to the terms B£ (a) and C A (b). 
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Fig. 4 (b) 



